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Abstract 

An analysis of orbit propagation models was performed by the Mission Operations element of the Sea-viewing 
Wide Field-of-view Sensor (SeaWiFS) Project, which has overall responsibility for the instrument scheduling. 
The orbit propagators selected for this analysis are widely available general perturbations models. The analysis 
includes both absolute accuracy determination and comparisons of different versions of the models. The results 
show that all of the models tested meet accuracy requirements for scheduling and data acquisition purposes. For 
internal Project use the SGP4 propagator, developed by the North American Air Defense (NORAD) Command, 
has been selected. This model includes atmospheric drag effects and, therefore, provides better accuracy. For 
High Resolution Picture Transmission (HRPT) ground stations, which have less stringent accuracy requirements, 
the publicly available Brouwer- Lyddane models are recommended. The SeaWiFS Project will make available 
portable source code for a version of this model developed by the Data Capture Facility (DCF). 


1. INTRODUCTION 

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) 
Project is a Code 970.2 activity at the National Aeronau- 
tics and Space Administration (NASA) Goddard Space 
Flight Center (GSFC). The primary responsibility of the 
Mission Operations element is to provide command sched- 
ules to maximize performance. Producing these command 
schedules requires the propagation of orbit positions over 
periods of days to weeks. It is of central importance to 
the successful performance of Mission Operations to select 
and use an orbit propagation model that will be accurate 
over these time scales. 

A secondary responsibility of the Mission Operations 
element is to provide pointing vectors and orbit models 
to worldwide ground stations to enable the acquisition of 
SeaWiFS data, which are directly broadcast in real time. 
These direct broadcast data are in High Resolution Pic- 
ture Transmission (HRPT) format. The time scales for 
these activities are typically much less than for command 
planning (usually less than three days). 

The purpose of this paper is to assess several widely 
available orbit propagation models for use by Mission Op- 
erations. Selection of orbit models for use by Mission Op- 
erations follows three steps: 1) collect requirements for 
performance, 2) assess availability of candidate models, 
and 3) analyze performance in relation to requirements. 
This analysis was a joint effort by the Mission Operations 
element and the SeaWiFS Data Capture Facility (DCF), 
which is responsible for supporting the HRPT ground sta- 
tions and has also implemented some of the tested orbit 
prediction models. 

2. BACKGROUND 

SeaWiFS is designed to make routine, global observa- 
tions of ocean color for a five-year mission lifetime. It is 
a follow-on sensor to the highly successful Coastal Zone 
Color Scanner (CZCS), which was carried on NIMBUS-7 


in operation from 1978 to 1986. The SeaWiFS instrument 
will be contained on the SeaStar spacecraft, on which it is 
the sole occupant. 

In a unique arrangement, the sensor, spacecraft, and 
launch vehicle (Pegasus) are being built and will be oper- 
ated by Orbital Sciences Corporation (OSC) of Chantilly, 
Virginia. NASA’s role will be to purchase the SeaWiFS 
data from OSC and provide command schedules to enable 
global coverage of the Earth and facilitate calibration of 
the sensor. Mission Operations serves as the component 
which determines these schedules and passes them to OSC 
for uploading. 

SeaStar will be placed in a 705 km, sun-synchronous, 
near-local-noon descending node orbit. These and other 
orbit characteristics are summarized in Table 1 . 


Table 1 . SeaStar spacecraft orbit parameters. 


Orbit Parameter 

Value 

Altitude 

705 km 

Eccentricity 

< 0.002 

Orbital Repeat Time 

16 days (233 orbits) 

Period 

98.9 minutes 

Inclination 

98.2° 

Equator Crossing Time 

Noon Local Time 

Node Type 

Descending 

Successive Equatorial 


Crossing Longitude 

-24.721° 


2.1 Mission Operations 

As stated earlier, the primary role of the Mission Op- 
erations element is to provide OSC with schedules of com- 
mand sequences to maximize the coverage and seientifi* 
usefulness of SeaWiFS data and ensure data acquisition. 
All of these commands are scheduled according to view- 
ing conditions or events which are a direct function of the 
spacecraft orbit position. The commands under the re- 
sponsibility of Mission Operat ions include: sensor eloctron- 
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ics on and off times, data down-link times at NASA’s Wal- 
lops Flight Facility (WFF), tilt change times, gain changes, 
calibration times, and broadcast times for research HRPT 
stations. 

2.2 Orbit Propagation Models 

There are basically two classes of orbit propagation 
models: general perturbations models (GPMs) and special 
perturbations models (SPMs). GPMs are reformulations 
of the equations of motion so that an analytical solution 
may be found. Typically, they include only the perturbing 
forces of a low r order Earth gravity field, although, some 
also may contain atmospheric drag. A common example is 
the Brouwer- Lyddane model. GPMs are fast and computa- 
tionally inexpensive at the cost of a reduction in accuracy. 

SPMs are characterized by the requirement of numeri- 
cal integration of the equations of motion. Such methods 
are computationally expensive and output is much slower 
than for GPMs. The advantage is that they can take 
into account more perturbing forces, such as a high-order 
Earth gravity field, multiple body gravitational forces, at- 
mospheric drag, solar radiation, tides, and others. Con- 
sequently, SPMs can be very accurate, but they can also 
be more unstable and are intimately sensitive to the initial 
conditions. 

2.3 Orbital Elements 

A spacecraft orbit can be parameterized in several ways. 
Three common representations are: mean element sets, os- 
culating elements, and orbit state vectors. Both mean and 
osculating element sets include the six classical Keplerian 
elements (semi-major axis, eccentricity, inclination, right 
ascension of the ascending node, argument of perigee and 
mean anomaly), but in fact they have quite different in- 
terpretations. Mean element sets represent an average of 
the orbit and are specifically designed for use with GPMs. 
Commonly distributed orbital elements sets by the U.S. 
Space Command (formerly the North American Air De- 
fense, or NORAD), Naval Space Surveillance (NAVSPA- 
SUR), and the National Oceanic and Atmospheric Admin- 
istration (NOAA) are all mean elements. Mean element 
sets are not used directly to compute the spacecraft po- 
sition, but are converted to osculating elements for this 
purpose. 

Osculating elements are an instantaneous Keplerian 
representation of the orbit; the osculating elements vary 
significantly over the orbit for low Earth spacecraft and 
are not readily propagated. Orbit state vectors are defined 
as the Cartesian position and velocity vectors in the geo- 
centric reference frame. Orbit state vectors are frequently 
used as the input to SPMs and are the standard output for 
all models. Orbit state vectors, or derived quantities such 
as latitude, longitude, and altitude, are the information 
required for scheduling. 


The SeaWiFS Project has two potential sources of or- 
bital parameters: distributed element sets such as those 
from NORAD or NAVSPASUR, or the Global Positioning 
System (GPS) data included in the spacecraft down link. 
While the GPS data is considered the primary source of 
orbit position information for navigation, at present there 
is no proven method for using this data in an orbit prop- 
agation model. Thus, distributed element sets are to be 
considered the first choice of orbit parameters. 

2.4 Requirements Analysis 

The most critical requirement for orbit propagation ac- 
curacy is to provide antenna pointing information for the 
down-link times at WFF. The Wallops antenna to be used 
for acquisition of SeaWiFS data is 9 in in diameter and 
autotracking with x-y tracking capability. The require- 
ments for acquisition of data by this antenna are 3 seconds 
along track, and 0.5° in azimuth, for a two-day propaga- 
tion. A lesser requirement is to provide predicted down- 
link times three weeks in advance for conflict analysis, for 
which the accuracy required is two minutes along track. 
The requirement for acquisition sets the limit for accuracy 
for the mission; none of the other commands have such a 
stringent accuracy requirement. 

Regarding HRPT ground stations, the SeaStar satellite 
L-band down-link has been designed to resemble all of the 
transmission parameters of the TIROS satellites. This was 
done in order to minimize the impact on small satellite 
ground stations around the world. The SeaWiFS Project, 
in turn, has adopted the TIROS data quality requirements, 
with a maximim bit error rate (BER) of 10“ 6 . In order for 
this BER to be met, OSC has suggested the system gain 
to total system noise temperature (G/T), a measure of a 
ground station’s quality performance, should be not less 
than 6.0. 

The antenna tracking requirements are given next. The 
L-band tranmission parameters are given in Table 2. 


Table 2. SeaStar L-band transmission parameters. 


Parameter 

Value 

Frequency 
Polarization 
Bandwidth 
Data Rate 

1207.5 MHz ±34.05 KHz 
Right-Hand Circular 
1.2MHz (— 3db) 
665.4 Kbps 


The pointing accuracies required for various recommended 
antenna dish sizes (including 1° of buffer) are shown in 
Table 3. 


Table 3. Pointing accuracies for different dish sizes. 


Dish Size [ft] 

Pointing Accuracy [°J 

8 

3.31 

6 

4.75 

5 

5.88 
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The SeaWiFS Project does not recommend using an an- 
tenna dish size of less than 5 ft since the BER of 10 -6 
is not guaranteed unless the system noise temperature is 
extremely low. 

These requirements are interpreted in terms of orbit 
prediction accuracies as follows. 

1. The 3-second timing requirement, for a typi- 
cal orbit velocity of 7.5 km s _1 , corresponds to 
22.5 km along track. 

2. The 0.5° azimuth error requirement at acqui- 
sion can also be evaluated in terms of along- 
track errors. At the planned SeaStar altitude, 
the station-to-satellite distance is approximately 
2,500 km at a typical acquisition of 5° above the 
horizon. 0.5° of azimuth error corresponds to 
21.8 km of error in the predicted satellite posi- 
tion perpendicular to the line of sight. 

3. For the HRPT stations, the 8-foot dish has the 
most stringent requirement at 3.31°; assuming 
that small stations do not have autotracking an- 
tennas, this requirement would have to be met 
by the orbit propagator throughout the contact. 
The orbit propagation errors would have max- 
imum effect when the satellite is directly over- 
head; at an altitude of 705 km, 3.31° corresponds 
to 40.7 km. 

Thus, the WFF acquisition timing and azimuth require- 
ments are most stringent and roughly equal in terms of 
orbit propagation accuracy. 

2,5 Approach 

The approach presented here is limited to assessing the 
feasibility of several widely available GPMs. The fast, com- 
putationally inexpensive analytical solutions of such mod- 
els, in conjunction with mean element sets, are appealing 
if the requirements can be met. 

3. METHODS 

3.1 Selection of Models for Evaluation 

The selection of orbit propagation models for evalua- 
tion was based on the following criteria: 

a) Availability of source code; this was necessary 
to allow the model to be configured as either a 
standalone program or a subroutine, to enable 
it to be run on any of several platforms and also 
to allow the format of the output to be closely 
tailored to the needs of Mission Operations; 

b) Use of known methods and models, with refer- 
ences to published derivations or specification; 
and 

c) Compatibility with readily available orbital ele- 
ment sets. 


For the last criterion, either NORAD two-line element 
sets or Navy elements from NAVSPASUR were considered 
as readily available for SeaWiFS. The so-called TBUS el- 
ements provided by NOAA for their satellites were not 
considered, because NOAA is not expected to provide this 
service for SeaWiFS. 

Two orbit propagation models were selected for eval- 
uation; each was available in two separate implementa- 
tions. The first was a standard Brouwer-Lyddane model, 
an analytic propagator which does not include atmospheric 
drag. The two implementations of this were as follows: 
a set of routines previously developed by one of the au- 
thors (Hoisington) for the DCF, using the original pub- 
lished work of Brouwer (1959) and Lyddane (1963), and 
also specifications for the Goddard Trajectory Determi- 
nation System (GTDS) (Cappellari et al. 1976); and the 
BRWLYD routine developed by NOAA (Kid well 1991). 
The latter model is presumably already widely used in the 
NOAA and Advanced Very High Resolution Radiometer 
(AVHRR) user communities. Brouwer-Lyddane propaga- 
tors use a mean element set as input; the derivation of this 
set from NORAD two-line elements is discussed below. 

The second model was the SGP4 low Earth orbit, model 
developed by NORAD Project Spacetrack. This is a mod- 
ified Brouwer propagator, which includes an atmospheric 
drag term, and is specifically designed to be compatible 
with NORAD element sets. This model was also imple- 
mented by Hoisington for the DCF from NORAD docu- 
mentation (Hoots and Roehrich 1980). After the evalua- 
tion was started, an updated version of the SGP4 source 
code was obtained directly from Project Spacetrack and 
was used for additional testing. 

As a side note, the SGP4 orbit model has restrictions 
placed on its use — it can only be used by U. S. government 
agencies and their contractors. 

3.2 Orbital Elements Used for Evaluation 

The evaluations were performed using NORAD two- 
line element sets for the following NOAA and Land Re- 
source Satellite (LANDSAT) spacecrafts: NOAA-12, which 
is in a polar orbit at a higher altitude than planned for Sea- 
Star (822 versus 705 km); and LANDSAT-4 and 5, which 
are in very similar orbits to that planned for SeaStar. Ele- 
ment sets were available at frequent intervals for all three 
satellites. 

The LANDSAT satellites perform regular orbit adjust- 
ment maneuvers, making extended evaluations impossi- 
ble, since no thrust model is available. NOAA-12 is com- 
pletely free flying, but at its higher altitude it is expected 
to experience a lower atmospheric drag effect. Thus, the 
NOAA-12 elements were used for extended studies, while 
the LANDSAT elements were used for short-term (up to 
10 days) analyses, as possible between orbit maneuvers, 
and also to validate the evaluation performed with the 
NOAA-12 elements. The actual element sets were obtained 
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from a bulletin board maintained by GSFC Code 513, and 
spanned the dates of July 1 through November 16, 1992. 

The use of the NORAD two-line element sets for the 
Brouwer-Lyddane models required a conversion since the 
NORAD elements are not explicitly mean element sets. 
Specifically, the mean motion provided with the NORAD 
elements sets must be converted to a semi-major axis. This 
conversion does not use the classical form, as described for 
Keplerian orbits in standard textbooks. The actual con- 
version was taken from the SGP4 source code and uses the 
mean motion, inclination, and eccentricity (a detailed de- 
scription is provided in Appendix A). Mean element sets 
derived in this manner were found to be entirely satisfac- 
tory for use in the Brouwer-Lyddane model. 

3,3 Constants Used in Models 

All of the models utilize constants to specify the gravi- 
tational field terms and Earth radius. In order to compare 
the various models and implementations, it was necessary 
to choose a consistent set of constants. The constants se- 
lected were taken from the system accepted by the Inter- 
national Astrophysical Union (IAU) in 1976 and published 
in the 1984 Astronomical Almanac (see Table 4). 


Table 4. Modeling constants for the Earth. 


Constant 

Value 

Radius, R e 

6,378.137 km 

Gravitational Constant, G e 

398,600.5 km 3 s~ 2 

Gravity Field Terms: 


J 2 

0.0010863 

J3 

-0.0000254 

J 4 

-0.0000161 


The Brouwer-Lyddane models also include a J 5 grav- 
ity field term, which was not specified in the IAU system 
of constants; both of the versions tested had this constant 
specified, but with different values. These constants were 
retained for the absolute accuracy analyses, but for com- 
parison purposes they were set to zero. 

3.4 Evaluation Methods 

The evaluation was performed in two parts: 1) determi- 
nation of the absolute propagation accuracy of each model, 
and 2) comparison of the propagations between the mod- 
els. The former was used to determine whether the models 
would meet the accuracy requirements for SeaWiFS Mis- 
sion Operations, while the latter was performed mainly to 
compare different implementations of the same model. 

The final output of each model is in the form of Carte- 
sian orbit state vectors (position and velocity) in the geo- 
centric inertial reference frame. The two parts involved 
a somewhat different approach to selecting the orbit vec- 
tors used for the comparisons, but in each case the actual 
vector comparisons were identical. 


Determining the absolute accuracy of a propagated or- 
bit requires the availability of a truth model. All of the 
orbit information for the subject spacecraft was obtained 
from the same source: the NORAD two-line elements sets. 
NORAD elements are intended to support scheduling re- 
quirements, as described earlier, not high-accuracy appli- 
cations (such as navigation). In fact they are believed to be 
degraded at the level of a few kilometers. This is still suffi- 
cient to determine whether the propagated orbit meets the 
SeaWJFS scheduling requirements, which are on the order 
of tens of kilometers. 

The use of NORAD elements at epoch as a truth model 
would be ideally validated by comparison with an exter- 
nal source of orbit data. While SeaWiFS Mission Opera- 
tions does not have access to independent data sources for 
the NOAA-12 and LANDS AT orbits used for this analysis, 
an evaluation of GPS data from the Extreme Ultraviolet 
Explorer (EUVE) satellite is currently being performed. 
The NORAD elements for EUVE are also available, but 
at longer intervals than for the NOAA-12 and LANDSAT 
elements. However, the GPS data, which spans the epoch 
time for a few of the NORAD element sets, demonstrates 
that the elements are accurate to within one kilometer at 
epoch. 

It was therefore assumed that the orbit vector deter- 
mined from each element set at its epoch time would be a 
reasonable truth model since no propagation was involved, 
just a conversion from mean elements to vectors. Thus, 
the approach in determining absolute propagation accu- 
racy was as follows. For each element set, an orbit state 
truth vector was computed at the epoch time. Then each 
element set was used to propagate the orbit to the epoch 
time of every other element set. This generated a number 
of orbit vectors representing various propagation times, all 
of which were compared with the truth vectors. 

The comparisons between models used a simpler ap- 
proach. For a given element set, orbit vectors were gen- 
erated at fixed intervals for a specified period using each 
model. In this type of evaluation, no particular set of vec- 
tors could be designated as truth, since all were generated 
by propagation; the point is to compare vectors from dif- 
ferent models which have been propagated for the same 
time since epoch. 

For both types of evaluations, the comparisons were 
performed as follows. The difference between two orbit 
position vectors was computed by subtracting each Carte- 
sian component (for the absolute accuracy evaluation the 
truth vector was subtracted from the propagated vector; 
for the comparisons the order of subtraction is purely ar- 
bitrary). The difference vectors were then converted to 
along track, cross track and radial (i.e., in the direction of 
the position vector) components. Determination of these 
components was performed by first determining the unit 
vectors in the directions of the orbit velocity, orbit normal 
(computed as the vector cross product of the position and 
velocity vectors), and orbit position. Specifically, for posi- 
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tion and velocity vectors P and V, and an orbit position 
difference vector D, 


and 


0 = PxV 


D c 




^rad 



(1) 

( 2 ) 

(3) 

(4) 


where £> at is the along-track position difference, D ci is the 
cross-track position difference, and D rac j is the radial posi- 
tion difference. 


4. RESULTS AND DISCUSSION 

The orbit propagation accuracy requirements for Sea- 
WiFS are stated for periods of 7 and 21 days. The results 
are shown for periods of 10 and 30 days, to demonstrate 
how well the requirements are met beyond the minimum 
times. Along-track propagation differences are expressed 
in both kilometers and seconds, since the primary accuracy 
requirements are stated in terms of timing errors. 

4.1 Absolute Accuracy Results 

This evaluation used 70 sets of NOAA-12 NORAD el- 
ements spanning the period July 1 through November 16, 
1992. The plots of absolute along-track propagation er- 
rors versus days since epoch for NOAA-12 are shown in 
Figs. 1-8. The results produced by SGP4, DCF/SGP4, 
DCF/Brouwer, and NOAA/Brouwer are shown for propa- 
gations of 10 and 30 days. 

The most significant factor in the absolute accuracy of 
the orbit propagations is the presence of a drag term in the 
SGP4 routines and the lack of such a term in the Brouwer- 
Lyddane routines. The drag model allows the mean along- 
track propagation error of the SGP4 routines to be near 
zero, even after 30 days, while the Brouwer propagations 
always show systematic (negative) errors, with what ap- 
pears to be a second order time dependence. 

All of the models met the stated timing accuracy re- 
quirements of 3 seconds after 2 days and 120 seconds after 
21 days. In fact, the 3-second requirement for acquisition 
was not exceeded for more than 7 days in all cases. The 
largest errors observed after 7 days were approximately 
10 km, or 1.5 seconds, for both of the SGP4 routines and 
15 km, or 2 seconds, for DCF/Brouwer. After 21 days, 
the largest errors were less than 15 seconds for all models. 
The plots for the SGP4 models showed a dispersion around 
roughly zero mean (Figs. 1-4), while the Brouwer models 
clearly showed the consistent degradation in accuracy from 
the lack of a drag term (Figs. 5-8). 


The azimuth requirement, equivalent to 21.8 km along 
track after 2 days, was also met by all of the models as 
shown by Figs. 1-8. The maximum errors after 7 days 
were approximately 10 km for the SGP4 models and 12 km 
for the Brouwer models. 

The most stringent HRPT station requirement, which 
is 40.7 km along track for an 8-foot dish, was easily met by 
all models well beyond 7 days. The Brouwer propagation 
errors were less than 40 km. for approximately 12 days, 
while the SGP4 models meet this requirement for at least 
two weeks. 

The other components of the propagation errors (cross 
track and radial) were consistently much smaller than the 
along-track errors, by as much as two orders of magni- 
tude. Sample plots of the cross-track and radial errors for 
DCF /SGP4 over 30 days are shown in Figs. 9 and 10. The 
maximum cross-track errors at 7 and 21 days were less than 
0.5 km and slightly more than 1 km, respectively, while the 
radial errors were less than 0.5 km and 1 km, respectively. 
Thus, the effects of orbit propagation on timing and an- 
tenna pointing can be assumed to be almost exclusively a 
result of the along-track errors. 

However, even with the SGP4 models, there are sub- 
stantial variations in the effectiveness of the drag model. 
This depends on the accuracy of the drag term included in 
the NORAD element sets. (This term is usually referred 
to as BSTAR in the documentation and code; no units are 
given.) Examination of the BSTAR term for the NOAA-12 
element sets showed large variations around the average, 
ranging from near zero to approximately double the aver- 
age value. The performance of the SGP4 propagation cor- 
relates very closely with these variations in the BSTAR term; 
the points on the SGP4 plots which show the largest neg- 
ative trend correspond to a BSTAR of approximately zero, 
and show essentially the same behavior as the Brouwer 
model (see Figs. 2 and 4). 

The cause of the large variations in BSTAR, which is 
determined by NORAD along with the other elements, is 
unknown. In practice, the extreme values of BSTAR can 
easily be found by a cursory examination of the elements 
and rejected; this was not done for this evaluation to avoid 
skewing the results in favor of the SGP4 model. Filtering 
of the elements sets in this manner would be expected to 
reduce the maximum SGP4 propagation errors to less than 
1 second at 7 days and less than 10 seconds at 21 days, 
respectively. 

As mentioned previously, the NOAA-12 orbit is higher 
than the altitude planned for SeaStar (approximately 822 
versus 705 km). The atmospheric drag at the higher al- 
titude would be expected to be lower, and therefore the 
propagation errors for NOAA-12 may not be considered 
truly representative for SeaStar. The LANDS AT-4 and 5 
orbits are much closer in altitude to SeaStar; however, due 
to the regular performance of orbit maintenance maneu- 
vers, it is not possible to perform 30-day evaluations using 
LANDSAT elements. An evaluation was performed using 
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Fig. 1. SGP4 aiong-track propagation errors for NOAA-12: 0-10 days. 
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Fig. 2. SGP4 aiong-track propagation errors for NOAA-12: 0-30 days. 
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Fig. 4. DCF/SGP4 along-track propagation errors for NOAA-12: 0-30 days 


SECONDS 





SQN003S 



KM 


F.S. Patt, C.M. Hoisington, W.W. Gregg, and P.L. Coronado 



o 2 4 6 8 10 

DAYS OF PROPAGATION 

Fig. 7. NOAA/Brouwer along-track propagation errors for NOAA-12: 0-10 days. 
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Fig. 8. NOAA/Brouwer along-track propagation errors for NOAA-12: 0-30 days. 
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Fig. 9. DCF/SGP4 cross-track propagation errors for NOAA-12: 0-30 days. 
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Fig. 10. DCF/SGP4 radial propagation errors for NOAA-12: 0-30 days 
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Fig. 11. SGP4 along-track propagation errors for LANDSAT-5: 0-10 days. 


LANDSAT-5 elements over 10 days to demonstrate the 
comparability with NOAA-12, and the results for SGP4, 
DCF/SGP4, and DCF/Brouwer are shown in Figs. 11- 
13. As these figures show, the 10-day performance for the 
SGP4 routines is very similar to that for NOAA-12, while 
the DCF/Brouwer errors are slightly larger (approximately 
2.5 km after 7 days versus 2 km for NOAA-12). Given that 
the NOAA-12 results met the 21-day requirement by a sub- 
stantial margin, it is expected the SeaWiFS propagations 
will meet this requirement as well. 

4.2 Comparisons 

Comparisons were performed of the different versions of 
each of the tw T o models: between the SGP4 and DCF /SGP4 
routines, and between the DCF and NOAA versions of 
the Brouwer-Lyddane model. Comparisons between SGP4 
and Brouwer-Lyddane models were not performed, since 
the absolute accuracy evaluations clearly showed signifi- 
cant differences in the outputs of the two models, almost 
entirely due to the presence of the drag term in SGP4. 

The comparison of the latest SGP4 routine with the 
DCF implemented version is of some interest. These rou- 
tines have a common heritage, i.e., NORAD Project Space- 
track. An examination of the code shows the latest SGP4 
implementation has been substantially updated from the 


version produced by DCF from the previous Project Space- 
track documentation. The question is whether the changes 
were in the logic or involved improvements to the model. 

The comparison was performed using a typical set of 
NOAA-12 NORAD elements. Orbit vectors were produced 
at regular (30-minute) intervals for 30 days. The differ- 
ences were computed as along-track and cross-track com- 
ponents. The results for the SGP4 models are shown in 
Figs. 14 and 15, which indicate that the propagation dif- 
ferences, while not completely negligible, are small com- 
pared to the overall propagation errors. The maximum 
along-track differences at 7 and 21 days are less than 0.2 
second and 0.5 second, respectively, while the maximum 
cross-track differences are less than 0.1 km and 0.2 km, re- 
spectively. The maximum cross- track differences show very 
linear behavior in magnitude with apparently zero mean, 
while the along-track differences show some periodic ef- 
fects with an overall linear trend and a positive mean. The 
along-track differences indicate that the mean of the ab- 
solute propagation accuracy distribution is closer to zero, 
possibly due to better drag modeling in the new SGP4. 
The cross-track differences appear to represent a small dif- 
ference in the precession rate of the orbit plane. 

The Brouwer model comparison (Figs. 16 and 17) show 
a very small linear trend in the along-track component 
(approximately 0.005 second after 21 days). The cross- 
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Fig. 12. DCF/SGP4 along-track propagation errors for LANDSAT-5: 0-10 days. 
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Fig. 14. Along-track propagation differences between SGP4 and DCF/SGP4. 
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Fig. 15. Cross-track propagation differences between SGP4 and DCF/SGP4. 
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Table Al. The format of the NORAD element set. 


A generic NORAD two-line element: 

Line 1: sssssUulyulduuAuyyddd . f f f f f f f f uu mm dluuuuuu m ®d2uuuuuubstar L jij L j0uu e l no 

Line 2: sssss uu ii . iiii u lll • llll u eeeeeee u www. wwuinim.iiiimmnuim.iiniiimiinnurrrrx 

Line 1 Variables 

Line 2 Variables 

sssss is the 5-digit spacecraft ID 

ii.iiii is the inclination (°) 

ly is the spacecraft launch year (not used) 

111.1111 is the right ascension of the 

Id is the spacecraft launch day (not used) 

ascending node (°) 

yy is the element epoch year (two digits) 

eeeeeee is the eccentricity (a decimal point 

ddd is the element epoch day-of-year 

is implied at the left) 

f fff f fff is the element epoch time of day 

www.wwww is the argument of perigee (°) 

mmdl is a mean motion time derivative (not used) 

mmm . mmmm is the mean anomaly (°) 

mmd2 is a mean motion time derivative (not used) 

nn.nnnnimnn is the mean motion (revs /day) 

bstar is the drag term (mantissa and exponent) 

rrrr is the orbit number (not used) 

elno is the element set number (not used) 

x is not used 

A typical NORAD two-line element set for NOAA-12: 

21263U U 91 U 32 UU A U 92183 . 30748338 uu . 00000257 uu 00000-0 uu 13333-3 u 0 u u3228 

21263 uu 98 . 6941 u 213 . 0068 u 0013736 u 146 . 

1404 u 214 . 0628 u 14 . 22063660 u 58773 


Note: The bstar term represents a value of 0.00013333 (0. 13333E-03) and the eccent ricity is 0.0013736. All other values 


are formatted explicitly. 


track component shows very similar behavior to the SGP4 
comparison, with maximum differences of approximately 
0.1km and 0.26 km at 7 days and 21 days, respectively. 
This also seems to indicate a small difference in the orbit 
plane precession rate. Given the overall performance of 
these models, the differences are too small to demonstrate 
an advantage for either model. 

5. RECOMMENDATIONS 

The evaluations have shown that any of the four prop- 
agators evaluated, in conjunction with NORAD elements, 
are capable of meeting the stated accuracy requirements 
for SeaWiFS scheduling purposes. The SGP4 propagators 
provide superior results, especially if rudimentary quality 
assurance of the element sets is performed. The major 
disadvantage of the SGP4 model is the restriction to gov- 
ernment agencies and contractors. The recommendations 
are, therefore, as follows: 

L The Mission Operations element should use the 
SGP4 model with NORAD elements to support 
scheduling. The performance of the new SGP4 
routine and the DCF/SGP4 version are very 
nearly equivalent for this purpose, but the new 
SGP4 appears to model drag better. 

2. The Brouwer-Lyddane model meets the needs 
of all non-government facilities which need Sea- 
WiFS orbit propagations, using NORAD ele- 
ments converted to mean elements as described 
in Appendix A. The requirements for this pur- 
pose are less stringent, as stated earlier. The 
SeaWiFS Project will provide the DCF Brouwer- 
Lyddane propagator to HRPT stations. 


Appendix A 


The NORAD two-line elements completely specify a space- 
craft's orbit; unfortunately they cannot be used directly with 
the Brouwer-Lyddane orbit model. Specifically, the NORAD 
model uses the mean motion, whereas the Brouwer model re- 
quires the semi-major axis. These two parameters are redun- 
dant for classical Keplerian orbits but the conversion is more 
complex for non-spherical gravitational fields (i.e., low Earth 
orbits) and mean element sets. 

The format of the NORAD element set is given in Table Al. 
The Brouwer models input the elements as a 6-element array, 
where the order of the elements is: seniimajor axis (km), ec- 
centricity, inclination, right ascension of the ascending node, 
argument of perigee, and mean anomaly. All but the semi- 
major axis can be copied directly from the NORAD element 
set. 


The conversion of the mean motion to the semi- major axis 
is performed as follows. First, the mean motion is converted 
from revolutions per day (n) to radians per minute (xno): 

xno = n mo (A1) 


The calculation of the semi-major axis uses the gravitational 
constant in units of fractional Earth radii ( Rl ‘ 5 ) per minute, 
and also the J 2 perturbation term. The Earth radius ( R e ) 
and gravitational constant (G e ) were defined earlier (Table 4) 
in units of kilometers and seconds, and J 2 (unitless) was also 
defined, as follows: 


R e = 6,378.137 km (A2) 

G e = 398,600.5 km 3 s -2 (A3) 

.72 = 0.00108263. (Al) 


The revised value of the gravitational constant used below 
(xke) is: 

xke = 60 \/ G c Re 3 (A5) 

= 0.0743668531 (A6) 
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where xke is in units of Earth radii 1 5 (R 1 5 ) per minute. The 
initial (classical) estimate of the semi-major axis (al) is: 


al 


/xke\ 2 / 3 

Vxno/ 


(A7) 


where al is in units of Earth radii. 

The perturbation corrections to the semi-major axis use the 
inclination (i), the eccentricity (e), and J 2 as follows: 


SPM Special Perturbations Model 

TBUS Not an acronym, but a NOAA orbit prediction 
message. 

TIROS Television Infrared Observation Satellite 
WFF Wallops Flight Facility 

Symbols 


temp = 
dell = 
aO = 


delO = 
aOdp — 


0.75 J 2 


3cos 2 (z) — 1 
(1 — e 2 ) 1 5 


temp 

Id 2 " 


al 


1 - dell 



dell j) + L^dellJ 

temp 
aO 2 
aO R € 


(1 -delO) 


(A8) 

(A9) 

(A10) 

(All) 

(A12) 


where aOdp is the mean semi-major axis in kilometers. This 
value is entered into the first location of the Brouwer element 
array. 

For the sample NOAA-12 element set listed in Table Al, 
n = 14.22063660 revolutions per day, i — 98.6941° and e = 
0.0013736. In this example, the equations above give a mean 
semi-major axis of 7,192.074 km. 


Glossary 

AVHRR Advanced Very High Resolution Radiometer 
BER Bit Error Rate 
CZCS Coastal Zone Color Scanner 
DCF Data Capture Facility 
EUVE Extreme Ultraviolet Explorer 

GPM General Perturbations Model 
GPS Global Positioning System 
GSFC Goddard Space Flight Center 
G/T System Gain/Total System Noise Temperature 
GTDS Goddard Trajectory Determination System 

HRPT High Resolution Picture Transmission 

IAU International Astrophysicai Union 

LANDSAT Land Resources Satellite 

NASA National Aeronautics and Space Administra- 
tion 

NAVSPASUR Naval Space Surface Surveillance 

NIMBUS Not an acronym, but a series of NASA experi- 
mental weather satellites containing a wide va- 
riety of atmosphere, ice, and ocean sensors. 
NOAA National Oceanic and Atmospheric Adminis- 
tration 

NORAD North American Air Defense (Command) 

OSC Orbital Sciences Corporation 
SeaWiFS Sea-viewing Wide Field-of-view Sensor 


aO Intermediate perturbation correction variable. 

al Orbital semi-major axis in units of Earth radii. 

aOdp The mean orbital semi-major axis in kilome- 
ters. 

D Orbit position difference vector. 

D a t Along-track position difference. 

D ct Cross-track position difference. 

D ra d Radial position difference. 
delO, dell Intermediate perturbation correction variables. 

e Orbital eccentricity. 

G e Gravitational constant of the Earth (398,600.5 
km 3 s -2 ). 

i Orbital inclination. 

J2 The J2 gravity field term (0.0010863). 

J 3 The J3 gravity field term (—0.0000254). 

J 4 The J4 gravity field term (-0.0000161). 

J5 The J 5 gravity field term. 

n Mean orbital motion in revolutions per day. 

O Orbit normal vector (P x V). 

P Orbit position vector. 

R e Mean radius of the Earth (6,378.137 km). 

temp Temporary perturbation correction variables. 

V Orbit velocity vector. 

xke Revised gravitational constant in units of Earth 
radii. 

xno Mean orbital motion in radians per minute. 
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